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Abstract 


During  the  period  of  January  1975  to  January  1082riMhniques  were  developed 
for  using  big^-frequency  (HF:  3-30  MHz)  radar  to  measure  ocean  currents  and  verti* 
cal  current  shears  in  the  upper  one  or  two  metres  of  the  ocean  surface.  The  measure¬ 
ment  of  current  velocity  with  an  HF  radar  depends  on  the  ability  of  the  radar  to 
precisely  measure  the  phase  velocity  and  direction  of  propagation  of  ocean  waves  whose 
wavelength  is  one-half  the  radar  wavel«xgth.  Ip  the  absence  of  a  current,  the  waves 
travel  at  a  speed  given  by  the  still-water  dispersion  rdation.  An  utiderlying  current  will 
modify  this  speed.  The  radar  measures  the  actual  phase  velocity  through  a  Doppler 
shift,  and  the  wavelength  of  the  ocean  wave  b  known  through  the  first-order  Bragg 
scattering  relation,  so  a  difference  between  observed  and  theoretical  still-water  phase 
velocity  can  be  calculated.  In  addition,  longer  ocean  waves  are  affected  by  currents  at 
deeper  depths  than  are  shorter  ocean  waves.  By  measuring  the  phase  velocity  at  several 
different  wavelengths,  it  b  possible  to  measure  a  vertical  current  shear  in  the  top  one 
or  two  metres  of  the  ocean  surface.  Thb  b  a  measurement  that  b  very  difficult  to  make 
by  any  other  means. 

A  portable  coherent,  pubed-Doppler  HF  radar  system  was  developed  -ai  Stanford 
and  used  in  several  experiments,  both  on  land  on  the  California  coast  and  on  board 
a  ship  during  the  JASIN  experiment.  Two  diffei'ent  antenna  systems  were  developed 
during  the  course  of  the  experiments.  The  land-based  experiments  used  a  balloon- 
supported  vertical  half-rhombic  transmnting  antenna  and  a  single  wideband  receiving 
loop.  The  shipborne  experiment  used  a\muiti-element  omnidirectional  transmitting 
antenna  with  a  steerable  phased  array  consbting  of  eight  wideband  loops  for  the 
receiver.  \ 

Land-based  eeperiments  conducted  at  Pescadcro,  California  demonstrated  that  a 
current  could  be  measured  by  an  HF  radar,  and  that  its  value  agreed  well  with  that 
measured  by  in-sita  drifting  spar  buoys.  In  addition,  there  was  evidence  of  a  vertical 
current  shear,  both  from  the  radar  measurements  and  from  the  buoy  measurements. 
The  equipment  and  data  processing  techiques  were  developed  during  the  Pescadcro 


experiments. 

The  JASIN  experiment  was  an  attempt  to  apply  these  techniques  to  the  measure¬ 
ment  of  surface  current  and  current  shear  in  the  open  ocean.  The  radar  system  was 
successfully  installed  on  board  a  ship  with  tolerable  interference  to  the  other  ship¬ 
board  activities.  The  steerable  antenna  was  quite  rugged  and  performed  as  expected. 
It  produced  antenna  patterns  consistent  with  the  physical  aperture  of  the  array.  The 
wind  velocity  during  the  JASIN  experiment  was  quite  low,  so  wind-  and  wave-generated 
currents  were  quite  small.  Nevertheless,  there  b  some  evidence  of  a  current  shear.  Its 
magnitude  b  small  and  near  the  resolution  limit  of  the  radar,  but  it  appears  to  be 
about  50%  higher  than  what  would  be  expected  from  the  wind  and  wave  conditions 
present  at  the  time  of  the  experiments. 
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Chapter  1 

Introduction 


This  report  covers  the  period  January  1975  to  January  1982.  The  aim  of  the 
research  described  here  was  to  develop  techniques  for  using  high-frequency  (HF :  3-30 
MHz)  radar  to  measure  ocean  currents  and  vertical  current  shears  in  the  upper  one  or 
two  metres  of  the  ocean  surface.  Some  of  the  research  was  conducted  using  land-based 
equipment  (at  a  site  on  the  California  coast  near  Pescadero),  and  some  was  conducted 
from  a  shipborne  platform  (the  R.V.  Atlantis  H  during  the  Joint  Air-Sea  Interaction 
Experiment  -  JASIN).  The  radar  equipment,  specially  developed  at  Stanford  for  this 
research,  was  the  same  in  both  cases,  but  the  antenna  systems  differed  widely  between 
the  two  experiments. 

The  measurement  of  current  velocity  with  an  HF  radar  depends  on  the  ability  of 
the  radar  to  precisely  measure  the  phase  velocity  and  direction  of  propagation  of  ocean 
waves  whose  wavelength  is  one-half  the  radar  wavelength.  This  is  first-order  scattering 
[Barrick,  1972a],  and  in  contrast  to  second-order  scattering  (Barriclc,  1971,  1972b], 
which  involves  two  or  more  ocean  wave  trains  often  travelling  in  different  directions, 
the  first-order  lines  are  usually  very  sharp,  often  only  a  few  millihertz  wide  at  a  radar 
carrier  frequency  of  several  megahertz.  The  Bragg  scattering  mechanism  has  been  well 
established  as  responsible  for  first-order  echoes  [Cromhie,  1955,  Barrick,  1972a|.  Thus  it 
is  possible  to  relate  the  phase  velocity  of  the  ocean  waves,  as  determined  by  the  Doppler 
shift  of  the  radar  echoes,  to  the  wavelength  of  the  ocean  waves,  as  determined  by  the 
wavelength  of  the  radar  signal.  The  difference  between  the  observed  phase  velocity  of 
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the  ocean  waves  and  its  value  predicted  by  the  deep-water  dispersion  relation  has  been 
attributed  to  surface  currents  [Stewart  and  Joy,  1074].  Furthermore,  this  difference 
often  is  a  function  of  ocean  (or  radar)  wavelength,  and  thb  is  attributed  to  a  vertical 
current  shear,  since  longer  wavelength  ocean  waves  are  influenced  by  currents  at  deeper 
depths  than  are  shorter  ocean  waves.  It  b  thb  wavelength  dependence  which  suggests 
that  an  HF  radar  can  be  used  to  measure  a  vertical  current  shear. 

The  first-order  scattering  mechanism  b  highly  selective  in  terms  of  the  direction  of 
wave  propagation:  only  those  ocean  waves  which  are  travelling  radially  toward  or  away 
from  the  radar  contribute  appreciable  energy  to  the  first-order  echo.  Ck>nsequently, 
ocean-wave  directional  resolution  b  obtained  by  radar  antenna  directivity.  Other  tech¬ 
niques,  such  as  those  used  by  the  CODAR,  [BarricJr,  Evans  and  Weber,  1077]  can 
provide  high  directional  resolution  with  smaller  antennas  than  those  used  here,  but 
they  were  not  employed  because  of  the  need  to  operate  over  a  wide  range  of  radar 
frequencies  in  order  to  measure  a  vertical  current  shear. 

Early  observations  were  carried  out  by  Stewart  and  Joy  [1074],  who  showed  that 
surface  currents  could  be  measured  by  an  HF  radar  system.  More  importantly,  they 
demonstrated  that  the  observed  current  was  a  function  of  radar  wavelength,  implying 
that  a  vertical  current  shear  could  be  measured  by  HF  radar.  To  determine  the  accuracy 
of  the  technique,  and  to  extend  these  early  observations,  we  conducted  a  series  of 
experiments  at  Pescadero,  California,  and  we  participated  in  the  JASIN  experiment. 
The  same  HF  radar  system  was  used  in  both  experiments,  but  with  different  antenna 
systems.  In  both  experiments,  the  Scripps  pitch-roll  buoy  [Stewart,  1077]  was  used 
to  measure  the  directional  ocean-wave  spectrum.  These  experiments  were  conducted 
jointly  by  Stanford  University  and  Scripps  Institution  of  Oceanography. 


Chapter  2 
Theory 


2.1  First-order  radar  scattering 

Crombie  [1055]  first  made  the  observation  that  radio  signals  scattered  by  the 
ocean  surface  consisted  of  a  pair  of  very  narrow  sidebands  dbplaced  above  and  below 
the  carrier  frequency  by  a  small  amount  (a  fraction  of  a  Hz  at  HF).  This  frequency 
corresponded  exactly  to  the  Doppler  shift  of  an  object  moving  at  the  phase  velocity  of 
an  ocean  wave  whose  wavelength  was  one-half  the  radio  wavelength.  He  deduced  that 
the  scattering  mechanism  was  Bragg  scattering,  in  which  the  contributions  from  each 
individual  wave  crest  all  add  up  coherently  to  produce  the  scattered  signal.  Hassebn&nn 
[1068]  showed  that  a  vertically-polarized  radio  wave  propagating  at  grazing  incidence, 
with  a  wavenumber  k*  and  frequency  (•;,■  is  scattered  into  two  waves  of  wavenumber  k* 
and  (jg  such  that 

k.  =  k.±k 

wp  =  w,  —  w,- =  iw  (1) 

where  uj)  is  the  Doppler  shift  of  the  scattered  radio  wave  and  k  and  u  are  the 
wavenumber  and  frequency  of  the  ocean  surface  wave.  For  the  case  of  backscatter 
(co-located  transmitter  and  receiver)  k«  =  — k,-,  so  that  2k«  =  i:k.  B&rrick  [1072b] 
examined  the  scattering  mechanism  in  more  detail  and  extended  the  geometry  to 
include  non-grazing  incidence. 
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2.  Theory 


The  ocean-wave  velocity  c  depends  on  an  underlying  current.  Let  the  still-water 
wave  velocity  be  ci  and  the  velocity  in  the  presence  of  a  current  be  Cp.  By  definition, 
Cp  =  (jjp/k  and  ci  =  oj/k.  Their  difference  v  is 

V  =  if'Cwp  —  w)  =  Ar“*  Aw  (2) 

where  Aw  is  the  difference  between  the  observed  Doppler  shift  and  its  expected  value. 
The  dependence  of  v  on  the  vertical  current  distribution  is  discussed  in  the  next  section. 
Stew&rt  and  Joy  [1974]  discuss  this  theory  in  more  detail. 

When  the  ocean  surface  roughness  is  calculated  to  second  order,  it  is  found  that 
the  Bragg  line  itself  has  sidebands.  B&nick  [1971]  showed  that  this  continuum  depends 
on  all  of  the  ocean  waves  present,  not  just  the  Bragg  component,  and  Hasselmann  [1971] 
suggested  that  this  can  be  used  to  measure  the  entire  ocean-wave  spectrum  using  a  single 
radio  wavelength.  However,  the  presence  of  the  second-order  sidebands  does  not  alter 
the  position  of  the  first-order  line,  so  as  long  as  its  position  can  be  accurately  deter¬ 
mined,  the  second-order  continuum  does  not  affect  the  current  measurement  capability 
of  the  radar  system.  Only  the  first-order  line  will  be  used  in  the  remainder  of  this 
report. 

2.2  Vertical  current  shear 

The  relationship  between  radar  wavelength  and  the  effective  depth  of  the  current 
measurement  was  first  examined  by  Stewart  and  Joy  and  has  recently  been  examined 
by  Ha  [1979]  as  part  of  a  Ph.D.  dissertation  at  Stanford  University.  The  difference  in 
velocity  between  a  wave  in  a  current  field  and  one  in  still  water  is 

roo 

U(k)  =  4kJ^  Uk{z)exp{—4kz)dz  (3) 

where  Uk(z)  is  the  component  of  the  horizontal  current  in  the  direction  of  the  radar 
beam  at  depth  z,  and  k  =  n/L  is  the  radio  wavenumber  for  a  radar  in  resonance  with 
odean  waves  of  length  L.  If  the  current  profile  is  logarithmic,  that  is 


(4) 
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where  u,  is  the  surface  drift,  u*  is  the  friction  velocity,  K  =  0.4  is  Von  Karman’s 
constant,  and  zq  is  the  reference  depth,  then  the  radar  probes  the  logarithmic  current 
at  a  depth  of  approximately  4.4%  of  the  ocean  wavelength  [Ha,  1070].  If  a  linear  profile 
is  assumed,  as  was  done  by  Stewart  and  Joy,  then  the  current  is  probed  at  a  depth  of 
approximately  7.8%  of  the  ocean  wavelength.  For  example,  in  the  JASIN  experiment 
the  ocean  wavelengths  were  31.2,  22.1,  11.2,  and  6.0  m.  For  this  case,  if  a  logarithmic 
profile  is  assumed,  the  radar  system  probes  the  current  at  depths  of  approximately  137, 
07,  40,  and  30  cm,  and  if  a  linear  profile  is  assumed,  the  depths  are  244,  173,  88,  and 
54  cm,  respectively.  Substituting  (4)  into  (3),  and  assuming  4Kzo  1, 

U{k)  —  +  \n{4k-tzo)  (5) 


where  n  =  1.781072.  Differentiating, 


dU{k)  ^  «• 
dlnF  K 


(e) 


where  F  b  the  radio  frequency  having  ocean  wavenumber  k  and  U{k)  b  the  current 
velocity  seen  by  thb  radio  frequency  [)Vu,  1075],  [Teague  et  a/.,  1077). 


2.3  Wind-induced  shear 


It  b  believed  that  the  wind  produces  a  vertical  current  shear,  because  the  wind 
blowing  over  the  sea  drags  the  surface  layers  along  in  the  direction  of  the  wind,  and  these 
surface  layers  further  influence  deeper  water.  In  order  to  estimate  the  magnitude  of  the 
wind-induced  shear,  it  b  assumed  that  the  velocity  profile  in  the  water  b  logarithmic 
mid  given  by  (4).  For  flows  past  a  free  boundary,  the  stress  r  b  given  by  r  —  ^pu%, 
where  p  b  the  fluid  density  and  ««  b  the  friction  velocity.  If  it  b  assumed  that  the 
stress  r  b  constant  across  the  air-sea  boundary,  then 


,  -  2  _  ^  _ 

•air  “  2^***®^***“^®**^  —  ^water 


2  _ 


(7) 
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2.  Theoiy 


If  the  stress  is  given  by  r  =  CppU^,  where  Cj}  =  I  ^  X  10~®  is  the  drag  coefficient  in 
the  air  and  U  is  the  air  velocity  at  10  m,  then 


(8) 


For  1/  =  10  ms  ^  thb  would  predict  a  current  shear  of  5  cm s  *  between  the  top  layer 
and  the  bottom  layer  probed  by  the  radar. 


2.4  Wave-induced  shear 

The  ocean  surface  waves  also  are  believed  to  produce  a  vertical  shear  through  the 
Stokes  drift.  The  Stokes  drift  arbes  because  a  wave  of  finite  amplitude,  accompanied 
by  irrotational  fluid  motion,  requires  that  the  vrater  particle  orbits  are  not  quite  closed; 
there  b  a  small  net  mass  transport  in  the  direction  of  wave  propagation  [Kinsman, 
1065,  p.  256].  Van  Dyke  [1082,  p.  110]  dbplays  some  striking  photographic  examples  of 
thb  effect.  In  order  to  estimate  thb  drift,  the  vertical  current  profile  b  predicted  using 
Stokes’  theory  of  wave-produced  currents  and  the  pitch-roll  buoy  measurements  of  the 
entire  set  of  ocean  surface  waves  present  during  the  radar  measurements.  The  vertical 
current  dbtribution  as  a  function  of  depth  z  b  estimated  from 

=  ^(w)w®exp(-2w*^/^)rfw  (0) 

where  g  b  the  gravitational  acceleration  and  F{u>)  b  the  surface  elevation  energy  density 
at  frequency  u  [Kenyon,  1060).  Then  for  each  of  the  waves  to  which  the  radar  system 
b  sensitive,  the  integrated  effect  of  thb  current  profile  on  the  wave  phase  velocity  is 
computed  by  using  the  current  distribution  (9)  in  (3).  In  performing  thb  calculation,  it 
b  assumed  that  the  Lagrangian  current  produced  by  a  wave  spectrum  advects  a  short 
wave,  thus  changing  its  phase  velocity.  Thb  is  not  precisely  true,  but  the  influence 
of  other  wave  interactions  b  small.  That  b,  the  influence  of  wave  interactions  is  to 
produce  a  current  which  modifies  the  phase  velocity  of  individual  wave  components; 


2.4.  Wave*ioduced  shear 


other  interactions  which  directly  change  the  wave  phase  velocity,  without  producing 
currents,  are  small.  The  shear  predicted  in  this  manner  is  on  the  order  of  l-6cms“‘  over 
the  depth  probed  by  the  radar  for  the  wave  conditions  encountered  in  the  experiments. 


Chapter  3 

HF  radar  system 


The  radar  system  used  in  the  experiments  was  developed  at  Stanford  specially  for 
this  program.  It  is  a  pulsed  coherent  Doppler  radar  which  can  transmit  and  receive 
pulse-interleaved  signals  simultaneously  on  four  different  frequencies  in  the  range  of 
4-30  MHz.  The  exact  frequencies  are  selectable  in  10  kHz  steps.  The  frequencies 
are  generated  by  a  set  of  phase-locked  loops.  The  same  signals  are  used  in  both 
the  transmitter  and  receiver,  thus  minimizing  the  amount  of  hardware  required  and 
ensuring  coherence  between  the  transmitter  and  receiver  internal  signals.  The  radar 
system  includes  a  digital  sampling  and  recording  system.  The  entire  radar,  exclusive 
of  the  power  amplifier,  occupies  12  inches  of  panel  space.  The  power  amplifier  occupies 
another  10  inches. 


3.1  TVansmitter 

The  transmitter  and  receiver  frequencies  are  generated  by  a  set  of  four  phase-locked 
loops.  These  each  generate  signals  in  the  range  of  34-60  MHz.  A  fifth  phase-locked  loop 
is  fixed  at  30  MHz.  The  transmit  signal  is  obtained  by  mixing  the  signal  from  one  of  the 
first  4  loops  with  the  30  MHz  signal,  producing  a  signal  in  the  range  of  4-30  MHz.  The 
signals  from  each  of  the  four  oscillators  arc  selected,  one  at  a  time  at  400  /ts  intervals, 
to  provide  the  interleaved  transmitter  pube  train.  The  pulse  widths  are  selectable  from 
10  to  200  /IS.  The  pulse  rise  and  fall  times  are  controlled  to  minimize  the  ocupied 


3.1.  Transmitter 


b.  Digital  system 


Figure  1 

Receiver  system  for  the  Stanford  four-frequency  radar. 

bandwidth.  The  signab  are  then  amplified  to  a  50  W  peak  power  level  by  a  commercial 
wideband  amplifier  to  form  the  final  transmitter  output.  All  of  the  phase-locked  loops 
are  referenced  to  a  stable  3  Mllz  crystal  oscillator.  The  use  of  the  double  loops  means 
that  the  original  8:1  frequency  range  can  be  covered  by  oscillators  which  need  cover 


10 


3.  HF  radar  system 


only  a  2:1  frequency  range. 

3.2  Receiver 

A  block  diagram  of  the  receiver  portion  of  the  radar  system  is  shown  in  Figure 
la.  The  receiver  consists  of  four  separate  amplifiers,  each  fixed-tuned  to  one  of  the 
transmitter  frequencies  and  followed  by  a  balanced  mixer  driven  by  one  of  the  same 
34-60  MHz  loops  used  to  drive  the  transmitter  (but  gated  on  for  a  full  400  /is).  The 
outputs  of  all  four  mixers  are  common.  This  generates  a  first  intermediate  frequency 
(IF)  of  30  MHz.  After  amplification  at  30  MHz,  the  signal  is  mixed  with  the  same  30 
MHz  transmitter  translation  loop  in  a  second  mixer  to  produced  a  second  IF  centered 
at  0  Hz.  The  use  of  in-phase  and  quadrature  local  oscillator  signals  at  the  second  mixer 
allows  the  separation  of  positive  and  negative  Doppler  shifts,  so  that  approaching  waves 
can  be  distinguished  from  receding  waves. 

3.3  Data  sampling  and  recording 

After  suitable  amplification  at  the  second  IF,  the  echoes  are  sampled  with  8-bit 
resolution  at  four  different  time  delays  (corresponding  to  four  different  radar  range  cells) 
after  the  transmit  pulse.  Both  the  time  delay  from  the  transmit  pulse  and  the  spacing 
between  samples  can  be  selected  in  10  /is  steps.  In  normal  operation,  the  samples  are 
taken  at  intervals  of  40  /is,  starting  50  /is  after  the  transmitted  pulse,  which  gives  a 
range  resolution  of  G  km.  A  number  of  such  complex  samples  are  added  together  in 
a  16-bit  accumulator,  both  to  improve  the  signal-to-nobe  ratio  and  to  allow  the  use 
of  a  low  peak  power  transmitter  by  using  a  high  pulse  repetition  rate.  The  number  of 
such  pulses  added  together  can  be  varied,  but  it  was  usually  set  so  that  the  effective 
sampling  rate  was  approximately  2  Hz.  This  allows  Doppler  frequencies  in  the  range  of 
—1  Hz  to  -Hi  Hz  to  be  resolved. 

The  digital  ])rocessor  is  shown  in  Figure  lb.  The  samples  are  accumulated  in  16 
separate  accumulators,  one  for  each  range  cell  at  each  frequency.  Each  one  is  32  bits 
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wide:  16  bits  for  the  real  part  and  16  bits  for  the  imaginary  part.  There  are  two  complete 
sets  of  buffers,  so  that  one  set  can  be  filled  while  the  other  set  is  being  recorded.  The 
resulting  16-bit  sums,  along  with  all  the  radar  system  paramaters,  are  recorded  in  digital 
form  as  a  pulse-code  modulation  (PCM)  waveform  on  a  1/4  inch  analog  instrumention 
tape  recorder.  For  the  shipbome  experiment,  the  receiving  antenna  was  composed  of  8 
elements  and  each  one  was  sampled  separately,  so  a  total  of  128  separate  complex  time 


Chapter  4 

The  Pescadero  Experiments 


The  Pescadero  experiments  were  designed  to  study  the  effects  of  vertical  current 
distribution,  determined  by  in-situ  measurements,  on  the  radar  signals.  Consequently, 
a  one-dimensional  geometry  was  adequate,  and  the  antenna  system  provided  a  narrow 
beam  for  all  four  frequencies  fixed  in  one  direction.  On  the  other  hand,  the  JASIN 
measurements  were  an  application  of  the  technique  to  obtain  oceanographic  data  in 
the  open  ocean,  and  some  degree  of  antenna  steering  capability  was  required.  As  a 
result,  the  antenna  systems  were  quite  different  for  the  two  experiments. 


4.1  Radar  installation 

On  several  occasions  between  May  1975  and  January  1078,  HF  radar  backscatter 
observations  of  the  ocean  surface  were  made  from  a  site  on  the  California  coast  near 
Pescadero.  A  map  of  the  site  of  the  experiment  is  shown  in  Figure  2.  The  site  was  located 
about  1  km  from  the  coast,  as  indicated  in  Figure  3.  For  these  experiments,  the  radar 
was  operated  on  frequencies  of  6.78,  13.38,  21.77  and  29.80  MHz.  These  frequencies  are 
in  Bragg  resonance  with  ocean  waves  of  length  22.1,  11.2,  6.9  and  5.0  m. 

The  transmitting  antenna  was  a  vertical  half-rhombic  antenna,  approximately 
250  m  long  and  45  m  high  at  the  apex,  supported  by  a  large  helium-filled  balloon.  This 
antenna  system  was  chosen  to  provide  a  moderately  directional,  wideband  transmitting 
antenna  at  a  reasonable  cost.  It  required  no  fixed  towers  to  install  and  maintain.  It 
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Site  of  the  Pescadero  experiment,  showing  the  locations  of  the  spar  buoys 
and  range  stations. 


did  have  the  disadvantage  that  the  balloon  had  to  be  filled  with  helium  before  each 
experiment,  so  it  required  about  an  hour  to  set  up  the  antenna.  Operation  of  the 
antenna  was  limited  to  wind  speeds  of  about  15  ms~'  at  the  radar. 


The  receiving  antenna  was  a  single  wide-band  air-core  loop  one  metre  in  diameter, 
manufactured  by  Hermes  Company.  The  loop  had  a  matched  preamplifier  installed 
right  at  the  loop.  The  combination  of  the  loop  antenna  and  matched  preamplifier 
resulted  in  an  antenna  which  had  an  effective  height  of  one  metre  over  the  entire  3-30 


n 


4.  The  Peseadero  Experiments 


Fipnre  S 

Propagation  paths  over  land.  This  is  a  rough  sketch  of  the  coastline;  small* 
scale  details  are  ignored. 


MHz  frequency  range.  Because  of  the  fixed  transmitting  antenna,  the  system  antenna 
pattern  was  fixed  about  a  nominal  center  direction  of  315” T.  The  3dB  full  width  of  the 
composite  radar  beam  varied  between  50”  and  12”,  as  shown  in  Figure  4.  The  measured 
beam  center  direction  varied  from  320”T  at  30  MHz  to  200”T  at  4  MHz  as  a  result  of 
the  propagation  path  over  land. 


Observations  were  usually  made  four  times  each  month  and  typically  consisted  of 
four  half-hour  runs,  each  with  a  different  radar  pulse  width.  Wind  speed  was  recorded 
locally  and  at  several  coastal  sites  operated  by  the  National  Weather  Service.  Several 
observations  which  included  both  the  Scripps  pitch-roll  buoy  [Stewart,  1077]  and  in-situ 
drifter  measurements  of  the  ocean  current  were  made  from  10-27  January  1078;  these 
observations  are  included  in  this  report.  The  locations  of  the  spar  buoys  and  the  radar 
range  cell  footprint  are  shown  in  Figure  2.  For  these  observations,  the  wind  speed  was 
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Figure  4 

Measured  3  dB  full  width  of  the  composite  Pescadero  antenna  vs.  radar 
operating  frequencies. 

also  recorded  on  the  boat  used  to  deploy  the  buoys. 

4.2  Spar  buoys 

The  drifters  were  four  lengths  of  spar  buoys  made  of  plastic  pipes  loaded  with 
suflicient  lead  shot  so  that  they  immersed  upright  in  water  to  predetermined  depths  of 
1,  3,  6,  and  12  m.  The  locations  of  the  buoys  were  determined  by  a  microwave  tracking 
system  which  had  a  range  resolution  of  a  few  metres.  The  boat  used  to  deploy  the  buoys 
carried  the  master  station  which  continuously  interrogated  two  transponders  located 
about  30  km  apart  on  the  coast.  The  position  of  each  buoy  was  measured  by  carefully 
maneuvering  the  boat  to  within  a  few  meters  of  the  buoy  and  then  noting  the  range 
to  each  of  the  two  transponders.  The  position  measurements  were  repeated  at  time 
intervals  of  an  hour  or  longer,  and  velocities  were  calculated  from  the  dilTcrenccs  in 
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Figure  5 

Smoothed  Doppler  spectra  at  four  radar  frequencies  and  four  range  cells. 
Data  were  collected  at  Pescadero  on  25  January  1978. 


position.  Because  the  buoys  typically  drifted  on  the  order  of  1  km,  the  uncertainty  of 
the  drift  velocity  was  less  1%  (Ha,  1079]. 

The  buoys  were  not  completely  immersed  in  the  water;  for  visibility,  a  portion 
of  the  buoy,  10-25%  of  the  buoy  length,  was  exposed  to  wind  drag  above  the  water 
surface.  Ha  calculates  that  the  effect  of  the  wind  drag  on  the  buoy  velocity  is  less  than 
the  statistical  uncertainty  in  the  buoy-inferred  drift.  He  also  demonstrates  that  a  buoy 
immersed  to  depth  d  will  drift  with  a  speed  equal  to  the  current  at  a  depth  of  O.Sd  and 
0.27(1  for  a  current  varying  linearly  and  logarithmically  with  depth,  respectively. 


To  obtain  an  indication  of  the  variation  of  ocean  current  within  a  single  radar  range 
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Figure  6 

Radar-inferred  depth-averaged  current  at  four  range  ceils  vs.  time  during 
25  January  1978.  Horizontal  lines  are  the  baselines  for  the  different  range 
cells. 

cell,  two  sets  of  buoys  separated  by  approximately  5  km  were  deployed  within  the  third 
range  cell  of  the  radar  system.  Because  only  one  boat  was  used,  the  measurements  of 
current  at  these  two  points  were  also  separated  by  30-60  minutes,  the  time  required  by 
the  boat  to  move  from  one  location  to  the  next  and  to  find  the  buoys. 
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Figure  7 

Comparison  of  radar-inferred  drift  and  spar-bupjr  drift  (cluster  2)  obtained 
on  19  January  1978.  The  horizontal  axis  is  logarithmic  in  depth,  and  a 
logarithmic  profile  is  assumed.  The  approaching  Bragg  line  was  used. 

4.3  Data  processing 

The  radar  echoes  at  four  different  frequencies  were  each  sampled  at  four  range 
cells  approximately  twice  a  second  for  a  total  period  of  2-4  hours  for  each  data  run. 
Fourier  transforms  were  computed  for  each  10-minute  data  segment.  The  first-order 
echoes  were  extracted  by  first  searching  for  the  peak  of  the  power  spectrum  in  a  narrow 
region  centered  on  the  theoretical  position  of  the  Bragg  line.  Then  the  null  between  the 
first-order  peak  and  the  second-order  continuum  on  either  side  of  the  Bragg  line  was 
found  by  searching  for  the  points  at  which  the  local  slope  of  the  smoothed  frequency 
spectrum  changed  sign.  The  nulls  thus  found  defined  the  limits  of  the  first-order  energy 
band.  The  centroid  of  the  first-order  spectrum  (the  first-order  moment)  was  found  by 
integrating  the  spectrum  over  this  range.  The  frequency  corresponding  to  the  centroid 
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Figurt  8 

Comparison  of  radar-inferred  drift  and  spar-buoy  drift  (cluster  1)  obtained 
on  25  January  1978.  The  horizontal  axis  is  logarithmic  in  depth,  and  a 
logarithmic  profile  is  assumed. 

of  the  spectrum  was  taken  as  the  measured  frequency  of  the  Bragg  line.  The  theoretical 
Bragg  frequency  was  subtracted  from  the  measured  frequency,  and  this  difference  was 
converted  to  a  current  velocity  by  t;  =  /X/2.  Details  of  the  data  processing  are  described 
by  Ha  [1078]  in  his  Ph.D.  thesis.  The  use  of  this  centroid  calculation  was  quite  effective 
because  the  antenna  beam  width  was  relatively  narrow,  as  noted  above. 

An  example  of  the  smoothed  radar  frequency  spectra  at  the  four  radar  frequencies 
and  four  range  cells  is  shown  in  Figure  5.  The  first-order  Bragg  lines  at  ir/a  used 
to  determine  the  current.  Note  that  the  nulls  between  the  first-order  lines  and  the 
second-order  continuum  are  well-defined.  They  usually  could  be  found  reliably  by  the 
algorithm  described  above.  The  altorighm  was  not  sensitive  to  the  land  echoes  centered 
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Figure  9 

Variance  in  drifts  of  spar  buoys  released  in  the  same  range  cell.  Cluster  1  was 
released  at  25.8  km  from  Pescadero,  and  cluster  2  was  released  at  20.5  km. 

Data  were  collected  on  25  January  1978. 

at  DC,  or  to  signals  beyond  a  narrow  region  centered  on  the  Bragg  line. 

4.4  Results 

Radar-inferred  currents  are  shown  in  Figure  6.  There  b  a  strong  time-dependent 
tidal  component,  but  in  addition,  there  b  variation  both  in  range  and  with  radar 
frequency.  Figures  7  and  8  show  the  comparbon  between  the  radar- inferred  current  and 
the  spar-buoy  drift  on  19  January  and  25  January  1978,  respectively.  For  these  figures, 
it  was  assumed  that  the  current  varies  logarithmically  with  depth,  an  assumption  which 
appears  to  be  justified  by  numerical  inversion  of  some  of  the  radar  data  |i/a,  1078], 
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Figure  10 

Logarithmic  current  shear  measured  by  the  four-frequency  radar.  The  horizon¬ 
tal  axis  is  logarithmic  in  depth,  and  a  logarithmic  profile  is  assumed  in 
the  inversion.  Vertical  bars  indicate  the  resolution  limited  by  a  finite-time 
Fourier  transform.  Data  were  obtained  at  Pescadero  at  approximately  14:15 
to  15:15  local  time;  range  6  to  13.5  km. 


The  radar  and  the  buoy  measurements  show  similar  changes  of  current  with  time 
which  are  apparently  related  to  the  tidal  cycle.  Also,  the  radar  and  buoy  measurements 
show  similar  decreases  of  current  with  depth.  On  25  January,  both  the  radar  and  buoy 
data  show  a  vertical  gradient  in  the  measured  current  which  increases  with  time.  This 
can  be  explained  in  part  by  a  rise  in  the  wind  speed  in  the  direction  toward  the  radar. 
However,  the  magnitude  of  the  vertical  shear  as  measured  by  the  radar  is  4-6  times 
the  value  predicted  by  (8)  from  the  wind  speed,  and  the  shear  measured  by  the  buoys 
is  also  higher  than  would  be  predicted  from  the  waves  alone. 


Figure  9  illustrates  the  variance  in  the  drifts  of  spar  buoys  released  in  the  same 


2i 


4.  The  Pescadero  Experiments 


Summary  of  Measured  and  Predicted  Shears  -  Pescadero 

Date,  1978 

Radar  meas. 

Total  est. 

Wind  est. 

Wave  est. 

19  Jan 

0.91 

4.81 

1.78 

3.03 

24  Jan 

4.69 

2.73 

1.55 

1.18 

24  Jan 

9.50 

2.73 

1.55 

1.18 

25  Jan 

1.07 

1.74 

1.40 

0.34 

25  Jan 

7.08 

1.74 

1.40 

0.34 

25  Jan 

9.80 

1.74 

1.40 

0.34 

27  Jan 

9.95 

4.95 

3.18 

1.77 

27  Jan 

13.31 

4.95 

3.18 

1.77 

Table  1 

Summary  of  Pescadero  radar  measurements  and  wind*  and  wave-predicted 
values  for  the  current  shear.  The  radar  values  represent  the  difference  in 
cm  s~*  between  the  current  velocities  obtained  from  radar  measurements  of 
5.0  m  and  22.1  m  ocean  wave  phase  velocities.  The  wind  and  wave  estimates 
of  the  shear  are  based  on  the  difference  in  current  velocities  predicted  for 
depths  of  21  cm  and  93  cm,  respectively.  These  are  the  depths  at  which  the 
5.0  m  and  22.1  m  waves  are  coupled  to  the  currents,  assuming  a  logarithmic 
current  profile. 

range  cell.  It  is  evident  that  the  two  sets  of  buoys,  separated  by  5  km  give  quite 
different  results.  In  fact,  by  comparing  Figures  8  and  9,  it  appears  that  thcft  is  better 
agreement  between  the  radar  measurements  and  either  set  of  buoy  measurements  than 
between  the  two  sets  of  buoy  measurements.  Thb  is  an  indication  of  the  horizontal 
variability  of  the  current  field.  Because  the  radar  illuminates  a  finite-sized  range  cell 
approximately  7.5  km  deep,  the  radar  measurement  represents  an  average  rather  than 
a  point  measurement,  so  the  two  methods  cannot  be  directly  compared. 

Figure  10  illustrates  the  current  shear  measured  by  the  radar  on  24  January,  and 
shows  that  a  logarithmic  variation  of  current  with  depth  provides  a  good  fit  to  the 
data.  The  slope  of  a  visual  fit  to  the  data  points  provides  an  estimate  of  1.6  cms~^  for 
«*,  a  value  about  twice  what  would  be  predicted  from  the  wind  speed  of  6  ms“*. 

Table  1  summarizes  the  vertical  current  shears  measured  by  the  HF  radar  system  in 
several  experiments  during  January  1978.  Also  listed  are  values  of  vertical  current  shear 
which  would  be  predicted  by  (8)  and  (3)  from  the  wind  and  wave  conditions  present  at 
the  times  of  the  experiment.  The  radar  measurements  showed  a  wide  variability,  but  in 
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Chapter  5 

The  JASIN  Experiment 


During  the  summer  of  1978,  the  Joint  Air-Sea  Interaction  (JASIN)  experiment  was 
conducted  in  the  North  Atlantic  Ocean  off  the  coast  of  Scotland.  The  overall  goal  of 
this  international  experiment  was  to  measure  the  momentum  transfer  between  the  wind 
and  the  waves.  Most  of  the  oceanographic  measurements  were  conventional,  depending 
on  tethered  buoys  supporting  current  meters,  thermistors,  and  anemometers;  pitch-roll 
buoys  for  measuring  wave  height  and  direction;  and  instruments  lowered  from  ships  for 
measuring  the  temperature,  salinity  and  density  of  the  deep  ocean  water.  In  addition, 
daily  observations  of  HF  radar  backscatter  were  made  throughout  the  experiment  to 
determine  if  a  vertical  current  shear  could  be  measured  from  a  platform  in  the  open 
ocean. 

5.1  Radar  installation 

The  four-frequency  radar  was  installed  on  the  R.  V.  Atlantis-II  and  measured 
velocities  relative  to  that  platform;  ship  drift  can  be  estimated  from  the  shipboard 
Loran-C  navigation  system.  The  ocean-wave  directional  spectrum  was  measured  simul¬ 
taneously  with  the  radar  observations  by  a  pitch-roll  buoy.  The  directional  spectrum 
of  the  ocean  waves  can  also  be  inferred  from  the  directional  distribution  of  the  echo 
energy  arriving  at  the  radar. 

The  transmiting  antenna  consisted  of  a  set  of  quarter-wave  vertical  radiators  and 
was  nominally  omnidirectional.  Figures  11  and  12  show  the  transmitting  antenna. 
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Figure  II 

View  of  the  base  of  the  JASIN  transmit  antennas.  There  were  four  quarter* 
wavelength  antennas,  all  fed  in  parallel  from  a  common  cable.  The  ship 
was  used  as  the  ground  plane.  The  elements  themselves  were  made  of  two- 
conductor  transmission  line  (with  both  leads  tied  together)  to  lower  the  Q 
of  the  antenna  so  that  the  50  ftaec  pulse  could  be  passed  without  severe 
distortion  caused  by  a  narrow-bandwidth  antenna. 
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Figure  IS 

View  of  port  side  of  Atlantis-II  at  the  dock  in  Glasgow.  The  receiving  array  is 
visible  along  the  side  of  the  ship  as  six  light-colored  loop  antennas  (seen  edge- 
on),  each  with  an  associated  preamplifier  immediately  above  the  antenna 
element.  The  transmitting  antenna  was  not  installed  at  this  time. 


System  directivity  was  provided  by  an  array  of  eight  broadband  loop  receiving  antennas, 
six  on  the  port  side  of  the  ship  and  two  on  the  starboard  side.  The  six  antennas  on 
the  port  side  arc  shown  in  Figure  13,  and  the  individual  elements  are  shown  in  more 
detail  in  Figure  14.  The  loop  elements,  approximately  1  m  high,  were  constructed  as 
half-loops  and  made  use  of  the  ship  hull  to  provide  an  image  half-loop.  A  wideband 
preamplifier  was  mounted  immediately  above  each  loop. 


The  loops  were  selected,  one  at  a  time,  by  an  electronic  switch  preceding  the  radar 
receiver.  Thus,  the  signals  from  the  individual  loop  elements  were  recorded  separately. 
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Figure  15 

Plots  of  the  combined  transmit-receive  pattern  at  13.38  MHz.  The  plots  are 
linear  in  power.  The  direction  of  the  bow  of  the  ship  is  indicated  by  the 
arrows.  The  beam  has  been  steered  in  ail  5  directions  from  the  same  data 
set  by  multiplying  the  signals  from  the  individual  elements  by  appropriate 
complex  weights.  The  patterns  were  measured  by  using  a  transponder  on 
the  Meteor  while  the  Atlantis-Il  wm  slowly  turned  in  a  circle. 
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measured  by  using  a  transponder  on  another  ship,  the  R.  V.  Meteor,  and  by  using  these 
data  the  combined  array  pattern  could  be  computed.  The  shape  of  the  array  pattern 
changed  as  the  beam  was  steered;  it  was  not  a  fixed  pattern  which  was  simply  rotated 
in  space.  Figure  15  shows  the  antenna  power  patterns  which  were  achieved  at  13  MIlz. 

The  radar  was  operated  on  frequencies  of  4.80,  6.78,  13.38,  and  21.77  MHz,  which 
are  in  Bragg  resonance  with  ocean  waves  of  length  31.2,  22.1, 11.2,  and  6.0  m.  Coherent 
processing  of  the  radar  signals  for  a  period  of  approximately  10  minutes  provides  a 
frequency  resolution  of  1.6  mHz,  from  which  the  Doppler  shift  of  the  radar  target  can 
be  measured  with  a  resolution  of  4.8  cms~*  for  the  longest  wavelength  and  1.2  cms~* 
for  the  shortest. 

5.2  Data  processing 

The  majority  of  the  data  procesing  was  concerned  with  estimating  the  current 
shear.  However,  in  a  related  study,  Ehbdalla  (1082]  examined  the  use  of  super>resolution 
spectral  estimation  techniques  to  estimate  the  ocean-wave  directional  spectrum. 

5.2.1  Estimating  the  current 

To  first  order,  the  radar  responds  only  to  those  ocean  waves  which  have  a  wavelength 
L  equal  to  one-half  the  radar  wavelength  X,  and  which  are  traveling  radially  toward 
or  away  from  the  radar.  The  phase  velocity  vp  of  these  Bragg-resonant  ocean  waves 
imparts  a  Doppler  shift  of  fp  —  2vp/\  to  the  radar  signals,  where  vp  =  ^gL/2ir, 
g  is  the  gravitational  acceleration  and  L  =  X/2.  Thb  Doppler  shift  is  independent 
of  the  direction  of  arrival  of  the  ocean  echoes,  because  the  radar  is  only  sensitive 
to  those  waves  which  are  travelling  exactally  radially  to  the  radar.  In  addition,  if 
there  is  a  current  present  with  a  component  vc  in  the  direction  of  the  radar  beam, 
then  if  the  radar  beam  is  sufficiently  narrow  there  is  an  additional  Doppler  Shift 
of  2vc7/X  due  to  this  current  component.  For  a  land-based  narrow-beam  radar,  the 
current  can  be  estimated  directly  from  the  measured  Doppler  shift  of  the  radar  echo. 
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Figure  16 

Geometry  for  a  radar  on  a  moving  platform.  The  ship  heading  is  v  and  the 
ship  is  drifting  in  a  direction  0  «iih  a  velocity  vr>  The  radar  is  receiving 
echoes  from  a  target  with  a  velocity  vc  from  a  direction  $.  The  antenna 
pattern  is  fixed  relative  to  the  ship,  and  the  echoes  are  seen  by  the  antenna 
at  a  relative  direction  d* 

However,  the  situation  is  more  complicated  if  the  radar  is  mounted  on  a  moving 
platform,  such  as  a  ship,  and  if  the  antenna  pattern  is  broad,  as  it  was  during  the 
JASIN  experiment.  A  typical  geometry  b  shown  in  Figure  16.  If  the  ship  b  moving  in  a 
direction  fJ  with  a  velocity  vt  with  respect  to  the  water,  there  b  an  additional  Doppler 
shift  of  {2vt/\)cos{0  — 0),  where  0  is  the  direction  from  which  the  radar  signal  arrives. 
(The  absolute  current  and  the  ship  drift  can  be  lumped  together;  it  b  only  the  relative 
motion  which  b  measured.)  In  the  absence  of  any  relative  motion  between  the  ship  and 
the  water,  the  frequency  spectrum  of  the  first-order  Bragg  echo  appears  as  two  very 
sharp  lines  at  ±/b;  the  ship  drift  spreads  these  into  two  bands  of  frequencies  (much  as 
a  synthetic  aperture  radar  does)  which  extend  ±2v2'/X  from  the  Bragg  center. 

The  shape  of  the  frequency  spectrum  between  these  limits  depends  on  the  direc¬ 
tional  distribution  of  ocean  energy,  the  shape  of  the  radar  antenna  patterns,  and  the 
Jacobian  of  the  transformation  between  angle-of- arrival  and  frequency.  For  a  given 
frequency  /  <  fma*  where  fmax  =  Svy/X,  the  angle  of  arrival  relative  to  the  ship  is 
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±cos  *(///ma*):  the  ±  arises  because  cos(tf)  is  an  even  function  of  0.  Defining 


S(0)  =  ocean- wave  directional  spectrum: 

energy  density  /  unit  angle  for  waves  arriving  from  0 
a(^)  =  antenna  power  pattern  in  direction  ^  from  the  bow  of  the  ship 
1/  =  heading  of  the  ship 
vx  =  magnitude  of  ship  drift  relative  to  water 
P  =  direction  of  ship  drift  relative  to  true  north 

/max  ~ 


then  one  would  expect  the  frequency  spectrum  of  the  radar  signal  to  be  given  by 


EAf)  =  E  - »)  ,  , 

where 

«.(/) = (-ir '  cos-'(///„„) + 


(11) 


with  the  beam  steered  in  the  direction. 

In  order  to  estimate  the  current  shear,  it  b  necessary  to  solve  equations  (10)  and 
(11)  for  vj  and  0  at  the  four  different  radar  wavelengths.  In  order  to  do  thb,  a  curve¬ 
fitting  procedure  was  used,  based  on  a  model  the  ocean- wave  directional  spectrum  as  a 
single-peaked  distribution  of  the  form  cos*((^  —  0o)/2).  The  values  for  the  ocean  model 
parameters  s  and  could  be  obtained  from  the  Scripps  pitch-roll  buoy  measurements, 
from  synthetic  aperture  HF  radar  data  taken  during  the  experiment,  or  left  as  free 
parameters.  Typically,  the  directional  dbtributions  during  JASIN  had  a  half-power 
beamwidth  of  30®  or  less,  much  narrower  than  were  found  by  Tyler  et  a/.  [1074]  at 
Wake  Island.  In  thb  case,  the  effective  direction  in  which  a  current  measurement  is 
made  depends  more  on  the  ocean  directional  dbtribution  than  on  the  antenna  pattern. 

The  antenna  beam  was  formed  in  five  different  directions  covering  a  00®  range,  and 
the  frequency  spectrum  of  the  radar  echo  was  computed  for  each  of  these  directions.  The 
frequency  estimate  E(f)  was  computed  for  each  of  these  directions  and  was  compared 
with  the  measured  spectrum  (the  shape  of  ay(^)  is  different  for  each  beam  direction 


j).  An  error  measure  over  all  five  directions  was  computed,  and  this  was  minimized  by 
varying  the  parameters  vt  and  ^  (and  in  some  cases  a  and  Bq  as  well)  until  a  minimum 
overall  error  was  achieved. 

Early  attempts  to  measure  the  vertical  current  shear  using  a  land-based  radar 
estimated  the  position  of  the  shifted  Bragg  line  by  noting  the  peak  of  the  frequency 
spectrum  [Steivart  and  Joy]  or  by  calculating  its  centroid  [Haj.  Such  procedures  in¬ 
herently  assume  that  the  radar  antenna  beamwidth  is  sufficiently  narrow  that  the  an¬ 
tenna  selects  only  echoes  (and  hence  current  components)  along  a  very  narrow  direction. 
Although  this  assumption  was  valid  for  the  Pescadero  experiments,  it  was  not  valid  for 
the  JASIN  experiment,  and  echoes  from  as  much  as  90®  can  contribute  to  the  observed 
frequency  spectrum.  In  this  case,  the  detailed  shape  of  the  frequency  spectrum  (and 
hence  its  centroid)  is  highly  dependent  on  the  shape  of  the  ocean  directional  spectrum 
and  the  antenna  pattern.  The  curve-fitting  procedure  described  above  was  originally 
implemented  with  a  gradient  search  over  the  parameters  of  the  models  for  both  the 
ocean  directional  spectrum  and  the  current,  but  the  data  were  too  noisy  for  this  ap¬ 
proach  to  be  stable.  Consequently,  the  model  parameters  for  the  ocean  directional 
spectrum  were  fixed  at  those  estimated  from  the  buoy  measurements,  and  the  current 
magnitude  and  direction  were  left  as  free  parameters  to  be  found  by  a  straightforward 
sequential  search.  Finally,  the  search  was  repeated  on  just  the  magnitude  of  the  current, 
with  its  direction  fixed  at  the  average  value  determined  in  the  last  step.  This  procedure 
was  applied  to  the  data  from  each  of  the  four  range  bins  at  each  of  the  four  radar 
frequencies. 

Once  the  apparent  current  was  obtained  for  each  radar  wavelength,  a  linear  regres¬ 
sion  analysis  was  applied  to  the  radar-inferred  current  as  a  function  of  the  logarithm 
of  the  ocean  wavlength.  The  shear  was  taken  as  the  differeuce  in  currents  given  by 
the  intercepts  of  this  line  at  the  the  minimum  and  maximum  depths  probed  by  the 
radar.  The  shear  value  was  retained  only  if  the  correlation  coefficient  of  this  fit  was 
greater  than  0.20,  and  only  if  at  least  3  of  the  4  radar  frequencies  were  used  (in  some 
cases,  individual  frequencies  were  omitted  because  of  strong  interference  identified  by 
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Figure  18 

Buoy,  ME,  ML,  and  MUSIC  spectrum  estimators.  5  Sept.  78,  6.78  MHz.  Tbe 
vertical  scale  is  linear  in  relative  power.  The  x>axis  represents  degrees  True. 
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Figure  19 

Buoy,  ME,  ML,  and  MUSIC  spectrum  estimators.  5  Sept.  78,  13.38  MHr. 
The  vertical  scale  is  linear  in  relative  power.  The  x-axis  represents  degrees 
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Figure  SO 

Buoy,  ME,  ML,  and  MUSIC  spectrum  estimators.  5  Sept.  78,  21.77  MHz. 
The  vertical  scale  is  linear  in  relative  power.  The  x^axis  represents  degrees 
True. 
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4.80  MHz 

Method 

Noi 

Nio 

Noi 

N20 

Nn 

Buoy 

0.666 

-0.335 

0.758 

0.432 

MUSIC 

0.507 

-0.528 

0.475 

0.526 

ML 

-0.141 

-0.393 

0.178 

0.822 

-0.024 

ME 

-0.571 

0.641 

0.377 

0.623 

-0.382 

6.78  MHz 

Method 

Noi 

Nio 

No2 

N20 

Nn 

Buoy 

0.583 

-0.586 

0.607 

0.648 

-0.237 

MUSIC 

0.602 

-0.650 

0.481 

0.510 

-0.395 

ML 

0.272 

-0.460 

0.454 

0.546 

-0.062 

ME 

-0.144 

0.330 

0.171 

0.820 

0.117 

13.38  MHz 

Method 

Not 

Nto 

No2 

N20 

Nn 

Buoy 

0.436 

-0.529 

0.813 

1.030 

0.039 

MUSIC 

0.800 

-0.010 

0.798 

0.202 

0.001 

ML 

0.705 

-0.110 

0.683 

0.316 

-0.055 

ME 

0.667 

0.238 

0.571 

0.429 

-0.007 

21.77  MHz 

Method 

A^oi 

Nw 

Noz 

N20 

Nn 

Buoy 

0.308 

-0.421 

0.568 

0.714 

0.068 

MUSIC 

0.601 

-0.386 

0.476 

0.524 

-0.164 

ML 

0.530 

0.472 

0.528 

-0.147 

ME 

0.933 

-0.150 

0.904 

0.095 

-0.104 

Table  2 

Comparison  of  moments,  5  Sept.  1078.  Wind  speed  12.3  ms~*  from  120°r. 


inspection  of  the  smoothed  radar  frequency  spectra).  A  linear  re^^cssion  analysis  was 
applied  individually  to  these  remaining  shear  values  as  a  function  of  the  shear  predicted 
by  (0)  for  the  waves  and  that  predicted  by  (8)  for  the  wind,  and  finally  a  multiple  linear 
regression  analysis  was  performed  against  both  the  wind-  and  wave-predicted  shears. 


Summary  of  Measured  and  Predicted  Shears  -  JASIN 

Date,  1978 

Radar  meas. 

Total  est. 

Wind  est. 

Wave  est. 

4  Aug 

6.8 

2.98 

1.34 

1.64 

5  Aug 


6  Aug 


7  Aug 


7  Aug 


8  Aug 


8  Aug 


10  Aug 


11  Aug 


11  Aug 


13  Aug 


13  Aug 


14  Aug 


24  Aug 


25  Aug 


25  Aug 


25  Aug 


28  Aug 


31  Aug 


31  Aug 


1  Sep 


1  Sep 


5  Sep 


5  Sep 


Summary  of  JASIN  radar  measurements  and  wind*  and  wave-predicted 
values  for  the  current  shear.  The  radar  values  represent  the  difference  in 
cm  s~*  between  the  current  velocities  obtained  from  radar  measurements  of 
6.9  m  and  31.2  m  ocean  wave  phase  velocities.  The  wind  and  wave  estimates 
of  the  shear  are  based  on  the  difference  in  current  velocities  predicted  for 
depths  of  30  cm  and  137  cm,  respectively.  These  are  the  depths  at  which  the 
0.9  m  and  31 .2  m  waves  are  coupled  to  the  currents,  assuming  a  logarithmic 
current  profile. 


5.2.2  Estimating  the  ocean-wave  directional  spectrum 


In  a  recent  Ph.D.  thesis,  Elabdalla  (1982]  examined  the  use  of  several  super¬ 
resolution  spectral  estimation  techniques  to  estimate  the  ocean-wave  directional  spectrum, 
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Figure  SI 

Comparison  of  ME,  ML,  and  MUSIC  spectrum  estimators.  This  figure  il> 
lustrates  the  improvement  using  the  modified  MUSIC  algorithm  as  com¬ 
pared  to  the  ML  and  ME  altorighms.  Data  were  collected  from  a  beacon 
transponder  located  at  —125".  The  vertical  axis  is  linear  in  relative  power. 
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and  he  applied  these  techniques  to  the  JASIN  HF  radar  data.  These  techniques,  nor¬ 
mally  used  to  estimate  a  frequency  spectrum  from  a  set  of  time  samples,  were  used  to 
estimate  the  ocean- wave  directional  spectrum  from  a  set  of  spatial  samples  at  the  posi¬ 
tions  of  the  receiving  loops  on  the  R.V.  Atl&Dtis-II.  He  compared  results  obtained  using 
the  Multiple  Signal  Classification  (MUSIC),  Maximum  Entropy  (ME),  and  Maximum 
Liklihood  (ML)  estimators  against  the  directional  spectrum  estimates  from  the  pitch- 
roll  buoy  and  synthetic  aperture  radar  (SAR)  processing  of  some  of  the  JASIN  data. 

In  order  to  demonstrate  the  directional  resolution  of  the  different  techniques,  the 
directional  spectrum  of  a  calibration  beacon  on  the  R.V.  Meteor  was  estimated.  Ideally 
thb  should  be  a  delta  function.  Figure  21  illustrates  the  results  of  this  test.  The  MUSIC 
algorithm  provided  an  unbiased  estimate  of  the  directional  spectrum,  while  the  ML 
and  ME  algorithms  provided  somewhat  biased  estimates.  All  three  altorighms  correctly 
identified  the  direction  of  the  beacon. 

The  algorithms  were  used  to  estimate  the  ocean-wave  directional  spectrum  under 
moderate  wind  conditions  (12  ms~^)  on  5  September  1078.  Figures  17-20  show  the 
results  of  these  algorithms  as  well  as  the  estimates  provided  by  the  pitch-roll  buoy. 
Again,  the  MUSIC  algorithm  gave  results  which  agreed  best  with  the  buoy  estimates. 
Table  2  compares  the  first  five  moments  Npf  for  the  various  algorithms.  Np^  is  defined 
by 

where  P(k,  ff)  b  the  ocean- wave  directional  spectrum  estimate  for  ocean  wavenumber 
k  and  direction  ff,  and  p  and  q  range  from  0  to  2.  Again,  the  moments  calculated  with 
the  MUSIC  algorithm  were  closest  to  the  buoy  moments.  Further  detaib  can  be  found 
in  lEkbdaUa,  1082]. 
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The  results  of  the  JASIN  experiment  are  summarized  in  Table  3  and  are  indicated 
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Wind-predicted  shear  (cm/s) 


Figure  it 

The  radar-measared  shear  vs.  the  shear  predicted  by  the  wind  alone.  The 
radar  shear  u  the  difference  in  the  deviations  of  the  phase  velocities  from 
their  still-waver  values  for  the  0.0  m  and  31.2  m  ocean  waves,  averaged  over 
all  range  bins.  The  squares  indicate  measurements  obtained  by  Ha  along  the 
California  coast  in  January  1078.  The  velocity  is  the  difference  in  current 
velocity  at  depths  of  30  cm  and  137  cm.  These  are  the  depths  at  which  the 
6.0  m  and  31.2  m  ocean  waves  are  effectively  coupled  to  the  currents  if  the 
current  distribution  is  lo^ithmic  with  depth,  ^m  eq.  (6),  these  can  be 
converted  to  by  dividing  by  3.78. 


5.3.  Results 


1  2  3  4  5  6  7 
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Figure  B3 

The  radar-measured  shear  vs.  the  shear  predicted  by  the  waves  alone.  The 
velocities  are  as  in  Fig.  22. 

graphically  in  Figures  22-24.  As  in  the  Pescadero  experiments,  the  estimate  of  the 
wind-  and  wave-generated  current  shear  was  made  with  the  assumption  of  a  logarithmic 
current  profile. 

The  wind  velocity  during  JASIN  was  generally  low,  reaching  10ms“*  on  only 


s  (cm/s) 
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Figure  S4 

The  radar>iceasured  shear  rs.  the  shear  predicted  by  the  sum  of  the  wind- 
and  wave-induced  components,  assuming  that  they  are  colinear.  The  velocities 
are  as  in  Fig.  22. 


two  or  three  occasions.  The  wave  heights  were  correspondingly  small.  Under  these 
conditions,  the  expected  shear  generated  by  the  wind  and  the  waves  are  nearly  identical, 
as  indicated  in  Table  3.  The  radar  measurements  indicated  vertical  shears  larger  than 
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predicted  by  either  the  wind  alone  (Figure  22)  or  the  waves  alone  (Figure  23),  but 
somewhat  less  than  would  be  expected  if  the  wind-  and  wave-induced  components  were 
assumed  colinear  and  added  together  (Figure  24).  For  this  data  set,  the  single-variable 
regression  on  the  waves  alone  resulted  in  a  slightly  higher  correlation  coefficient  than 
the  single-variable  regression  on  the  wind  alone  (0.71  vs.  0.70),  although  the  difference 
is  probably  not  significant.  The  two-variable  regression  on  both  the  wave-  and  wind- 
induced  shears  was  slightly  better,  resulting  in  a  correlation  coefficient  of  0.72,  with 
the  regression  equation 


—  (1.05  i  0.80)tlittt,at,e«  H*  (0.59  i  0.71)U;^u„‘,)f(  +  (1.54  i  1.18)  (f3) 

where  u*radar  is  calculated  from  the  radar-measured  shear  through  (6),  ^ 

timated  from  (8),  and  u*wavet  is  estimated  from  (0),  (3)and  (6).  The  standard  errors 
are  based  on  a  10%  risk  level. 
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Summary  and  Conclusions 


The  Pescadero  experiments  demonstrated  that  a  current  could  be  measured  by 
an  HF  radar,  and  that  its  value  agreed  reasonably  well  with  that  measured  by  in- 
situ  drifting  spar  buoys.  In  addition,  there  was  evidence  of  a  vertical  current  shear, 
both  from  the  radar  measurements  and  from  the  buoy  measurements.  The  design  of 
the  equipment  was  proven  and  the  data  processing  techniques  were  developed  in  the 
Pescadero  experiments. 

The  JASIN  experiment  was  an  attempt  to  apply  the  techniques  developed  in  the 
Pescadero  experiments  to  the  measurement  of  surface  current  and  current  shear  in  the 
open  ocean.  The  radar  system  offered  the  possiblity  of  making  a  measurement  that 
no  other  instrument  could  make.  The  radar  was  successfully  installed  and  operated 
on  board  a  ship  with  tolerable  interference  to  the  normal  shipboard  activities.  The 
steerable  antenna  proved  to  be  quite  rugged  and  successful.  The  antenna  patterns  which 
were  achieved  were  what  would  be  expected  from  the  physical  aperture  of  the  receiving 
array. 

The  v/ind  velocity  during  all  of  the  JASIN  experiment,  except  during  the  period 
when  the  ship  was  in  port,  was  much  lower  than  had  been  expected,  so  wind-  and  wave¬ 
generated  currents  were  quite  small.  Nevertheless,  there  is  some  evidence  of  a  current 
shear.  Its  magnitude  is  small,  however,  and  it  is  near  the  resolution  limit  of  the  radar 
system.  There  were  many  steps  involved  in  the  data  processing,  including  the  Fourier 
processing  of  the  radar  echo  signal,  the  formation  of  a  narrow  antenna  beam  in  several 
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different  directions  at  each  radar  frequency,  and  the  estimation  of  the  Doppler  shift  of 
the  radar  echo  by  the  estimation  of  four  parameters  of  a  model  for  the  ocean  directional 
spectrum  (center  and  width)  and  relative  motion  between  ship  and  water  (magnitude 
and  direction).  Consequently,  it  is  very  difiScult  to  compute  the  overall  precision  of 
the  radar  measurements.  The  internal  consistency  of  the  radar  measurements  from  one 
range  cell  to  the  next  suggests  a  figure  of  about  5cms“*.  Although  the  two-variable 
regresssion  indicates  that  the  wave-induced  component  of  shear  may  have  a  slightly 
higher  correlation  with  the  radar  measurements  than  the  wind-induced  component  (and 
thb  was  also  the  case  when  individual  range  bins  were  processed  separately),  this  may 
be  the  result  of  noisy  data.  However,  the  radar  measurements  indicate  a  shear  higher 
than  would  be  expected  from  either  component  alone.  Ha  also  noted  that  the  shear  he 
measured  at  Fescadcro  gave  a  value  about  twice  that  predicted  from  the  wind  alone 
(Figure  10). 
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